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En los sistemas de riego es necesario definir la línea de 
energía total del flujo de agua para evitar variaciones 
de presión y gasto en los puntos de entrega y control. En 
estos puntos se instalan equipos especiales, entre los que 
se encuentran los aforadores Venturi. Estos dispositivos 
han sido poco estudiados en relación con tamaños, 
formas, materiales o condiciones de funcionamiento, y las 
recomendaciones de operación provienen de características 
obtenidas de modo experimental. Así, modelar su 
funcionamiento a través de simulación con dinámica de 
Muidos computacional (DFC) resultaría práctico, económico 
y confiable. El objetivo del trabajo consistió en validar la 
simulación de funcionamiento con DEC de un dispositivo 
Venturi con base en observaciones de laboratorio. La 
comparación de los resultados experimentales y mediante 
DEC fueron altamente satisfactorios para los valores de 
velocidad media (153 m 5) gasto (0.027 m's') y presión en 
el manómetro diferencial (15 em de Hg). Una vez validado 
el modelo, se simularon ocho escenarios de operación, con 
variación de gasto desde 0.005 mv s' hasta 0.040 mm! s”, lo 
cual simpliicó el modelo de la relación gasto-pérdida de 
carga del Venturi, con una ecuación cuadrática y así eliminar 
la incertidumbre del coeficiente de descarga requerido en 
los análisis experimentales. Los valores característicos del 
dispositivo fueron un requerimiento de carga hidráulica 
mínima en la contracción del Venturi de 04 kg em”, y una 
pérdida de carga total en el Venturi de 0.075 kg cm* para 
Op. de 0.040 m? s!. Por último se recomienda diseñar 
dispositivos Venturi con DEC y con ello obtener información 
confiable para su instalación en sistemas de riego. 


Palabras clave: dinámica de fuidos computacional (DEC), 
aforador Venturi, sistemas de riego. 
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Irrigation systems. require defining the total energy line of the 
ont flow in order to prevent variations ón pressure and flow at 
delivery and control points. Special quipaner are installed at these 
points, including Venturi guuges. These devices ae not been 
ruidely studied in terms of sizes, slupes, materials or functioning. 
conditions, and the operating recommendations are derived from 
experimentally obtained. characteritic. Tus, modeling “tetr 
Fnctoning rot computationalfuid dypamics (CED) simudations. 
“would be practica, economica! and reliable. The objective of Mis 
ork 098 to validate the CFD simulation of the functioning of a 
Venturi decice based on Taboratory observations. The comparison 
of the experimental results and those obtained with the CED were 
Zig atifactory for the mean velocity (1.53 1 5), flow (0.027 
e 5) and pressure in the diferentil manometer (15 cm de Hg) 
Aer the model ws validated, operating cenarios vere sómulated 
vit variations in flow from 0.005 me 5* to 0.040 mé s* which 
simplified the modeling of the head loss fcw rate relationship of 
He Venturi, sit a quadratic equation, thereby climinating He 
vncertainty related to the discharge cocfícient required by te 
experimental analysis, The characteistic values of He device were 
A pninimcs hydraudic head of 0. kg cos at Ho contraction of the 
Venturi ana total pressure lossof0.075 kg cn” n the Venturi ora 
O, 00.040 0 5”. Lastl, the design of Venturi devices using CED. 
ds recommended to obtain reliable information for is installation in 
irvigaion systems 
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Introducción 


En el diseño y operación de estructuras hi- 
dráulicas integradas en plantas de bombeo, 
plantas potabilizadoras, centrales de generación 
hidroeléctrica, micro centrales hidroeléctricas y 
sistemas de riego, entre otras, es necesario co- 
nocer con la mayor precisión posible la línea de 
energía total y, en su caso, regularla, para evitar 
variaciones no deseadas de presión y gasto en 
los puntos de entrega y control. En efecto, estos 
puntos son hidráulicamente idóneos para insta- 
lar piezas especiales de medición de ujo, entre 
las que se encuentran los aforadores tipo Ventu- 
ri. Sin embargo, estos dispositivos han sido poco 
estudiados en relación con tamaños, formas, 
condiciones de funcionamiento y materiales, 
en particular debido a la elevada inversión que 
representa la experimentación en laboratorio de 
cada una de sus numerosas alternativas. 


Para el diseño de un aforador instalado en 
una tubería se utilizan conceptos y modelos de 
la mecánica de fluidos (Levy, 1957; White, 1994), 
lo que requiere definir, en primer término, el 
gasto por conducir (SARH, 1973). Post 
te se selecciona el tamaño idóneo de lon; 
y se definen las condiciones de alojamiento y 
operación del dispositivo. Bajo este proceso se 
revisan las dimensiones del aforador, conside- 


rando las corrientes cinéticas, así como las ca- 
ísticas del tipo de material y la velocidad 
interna del agua. 

La velocidad de un fluido a través de una 
tubería se puede determinar cuando se restrin- 
ge el flujo y se mide la reducción en presión 
debida al aumento de velocidad en el sitio de 
restricción. Los medidores de flujo que se basan 
en este principio se conocen como flujómetros 
de presión diferencial. Algunos de estos dispo- 
tivos de medición son los Venturi, las placas 
de orificio y las toberas, entre otros, y consisten 
básicamente de un elemento que estrangula al 
flujo y crea un cambio en la carga piezométrica, 
que se traduce en una pérdida de energía 

Chow (1959) señala que para definir las 
recomendaciones de operación de un equipo 
hidráulico, en particular en el estudio de piezas 


ractel 
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especiales como los aforadores Venturi (White, 
1994), es necesario determinar las pérdidas de 
carga entre los puntos de interés, Para estos 


cálculos se utiliza la ecuación de la conservación 
de energía (ecuación (1) 


(1) 


Donde v = velocidad del agua (m s”); 
3 = aceleración de la gravedad (m s?); z = carga 


de posición (m); Y, 


energía entre las secciones 1 y 2 (m); 2 = carga 


de presión en el punto a (m); E 
velocidad en el punto 1 (m), y y EN específico 
del fluido (kg m7). 

La ecuación es dedu 
y se utiliza en la mayoría de aplicaciones para 
flujo con variables promedio, como la veloci 


pérdidas totales de 


carga de 


dad media (m 5") en la sección considerada, La 
aplicación de este modelo se encuentra en todo 
tipo de estructuras hidráulicas utilizadas en el 
campo de la ingeniería (SARH, 1973). 

El aforador Venturi puede dividirse en 
cuatro zonas para su análisis (figura 1): a) zona 
de entrada: tiene como función favorecer una dis- 
tribución uniforme de las líneas de flujo dentro 
de la unidad y con ello la velocidad de llegada 
a la zona de contracción; b) zona de contracción: 


parte de la estructura en la cual se realiza el 
proceso de reducción de dí 


al diámetro menor por acción del estrangula- 
miento; c) zona de expansión: conformada por 
un tramo de transición del diámetro menor al 
diámetro mayor igual al de entrada; d) zona de 
salida: constituida por un tramo recto sin cambio 
de diámetro que conduce el lujo hacia la tubería 
en uso. 

Los componentes del aforador Venturi se 
diseñan para diferentes tamaños, formas, condi- 
ciones de funcionamiento y tipo de material. La 
zona de contracción-expansión puede ser curva 
o recta, según el diseño propuesto; así, para ca- 
da equipo se propone el ángulo determinante 
de la longitud de la ampliación, siendo justo la 
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saria para normalizar el flujo al final del 
dispositivo. Una de las relaciones más críticas 
en los Venturis es el grado de estrangulamiento 
(m1), que relaciona las áreas de la sección contraí- 
da (A,) con la sección de entrada (A); es decir, 
m =A,/A,. Para la medición de la pérdida de 
energía por el estrangulamiento, se recomien- 
da realizarla con base en dos secciones (figura 
1): una justo al final de D, y la segunda en la 
contracción D,, justo después de la terminación 
de contracción e inicio de la ampliación (ASME, 
1983). Para esta medición de la caída de la línea 
de energía entre las dos secciones señaladas 
existe una metodología de cálculo de acuerdo 
con el equipo disponible en laboratorio, sin 
embargo, el medio más recomendable es el 
manómetro diferencial de mercurio instalado 
en los puntos indicados. 

Por otro lado, la ASME (1983) presenta un 
modelo matemático deducido a partir de la 
ecuación (1) para la determinación del gasto 
en función de un coeficiente de descarga, con 
el conocimiento de la línea de distribución 
de presiones o línea piezométrica a partir 
del centro del tubo. Además, dado que la 
pérdida de energía entre las secciones de en- 
trada y de contracción está representado por 
h=(2,+0,/28)-(2,+03/28) y si lt se sustituye 
en la ecuación de conservación de la masa 
A, = A,v,/u, se obtiene la ecuación (2): 


Q=C,A 


Donde el coeficiente C, es el conocido coefi- 
ciente de descarga y se determina a partir de 
pruebas experimentales de laboratorio. Sme- 
tana (1957), reportado por Sotélo-Ávila (1979), 
presenta dos tablas para la determinación del 
coeficiente C; una en función de la relación de 
estrangulamiento m = A,/ A, y otra en función 
del número de Reynolds presente en la sección 
del estrangulamiento. 

Ranald (1978) y la ASME (1983) reportan 
diferentes tablas y figuras auxiliares en el 
cálculo de aforos de flujo en tuberías bajo el 


empleo de dispositivos tipo Venturi. Aquí cabe 
destacar que ambos presentan valores muy 
similares de coeficientes de descarga. Por su 
parte, la USDI (1979) y la 1OS (1991) reportan la 
estandarización de procedimientos que deben 
satisfacerse para la correcta medición de flujos 
con aforadores tipo Vent 

Tradicionalmente, para el cálculo de los pa- 
rámetros característicos del flujo, y debido a la 
complejidad de los fenómenos hidráulicos, era 
necesario discretizar el sistema de ecuaciones 
dejando fijas ciertas características del fuido, 
para después aplicar las ecuaciones generales de 
balance de masa y energía a un volumen finito 
del flujo y, en particular, plantear la ecuación que 
describiese su movimiento. Sin embargo, en la 
mayoría de los casos, esta descripción resultaba 
puntual y no permitía la representación espacial 
del sistema. Este inconveniente ha sido resuel- 
to con la dinámica de fluidos computacional 
(DFC) (Anderson, 1995). La herramienta DFC 
permite la visualización espacial del fenómeno 
en estudio, al obtener una solución numérica a 
los problemas de movimiento de fluidos, con 
un alto grado de precisión una vez validado el 
modelo planteado con datos experimentales, lo 
cual ha permitido una mejor comprensión de 
los fenómenos hidráulicos y de la mecánica de 
fluidos en general. 

En realidad, tanto la DFC y el análisis experi- 
mental son métodos complementarios, resaltan- 
do que el empleo de la DFC reduce la cantidad 
necesaria de análisis experimental en laborato- 
sio. Por ejemplo, en ingeniería hidráulica, se ha 
hecho investigación en el análisis de válvulas 
(Davies é: Stewart, 2002), bombas (Zheng, 2000) 
o enla aplicación de canales abiertos (Wu, Rodi, 
se Wenka, 2000). En agricultura protegida, el 
análisis con DFC ha sido empleado para mode- 
lar el clima al interior del invernadero y con ello 
plantear estrategias para su gestión mediante 
sistemas auxiliares de refrigeración o calefacción 
(Flores-Velázquez, Mejia, Rojano, € Montero, 
2011). En sistemas de riego presurizado es 
posible encontrar aplicaciones para el análisis 
de emisores de riego (Palau, Arviza, Se Frankel, 
2004) y dispositivos de inyección de fertilizantes 
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(Manzano 6 Palau, 2005). 

Diversos autores han simulado el compor- 
tamiento hidráulico de vertedores y aforadores 
usando técnicas DFC, después de validarlos 
con datos experimentales. Gharahjeh, Aydin y 
Altan-Sakarya (2015) simularon el flujo de agua 
para un vertedor de cresta delgada, mientras 
que Naghavi, Esmaili, Yazdi y Vahid (2011) 
simularon un vertedor circular; Temeepattana- 
pongsa, Merkley, Barfuss y Smith (2014) vali- 
daron un modelo DFC y ajustaron ecuaciones 
algebraicas genéricas para una diversidad de 
aforadores sin garganta “Cutthroat flume”, a 
partir de datos simulados de carga-gasto bajo 
condiciones de descarga libre y ahogada. 

El planteamiento de la DFC se fundamenta 
en la solución de las ecuaciones que definen el 
movimiento de un fluido, combinando dinámica 
de fluidos, programación y métodos numéricos. 
En general, el movimiento de un fluido está 
basado principalmente en procesos físicos, que 
pueden ser planteados en términos matemáticos 
como una serie de ecuaciones en derivadas par- 
ciales, mismas que representan numéricamente 
las componentes de un flujo. Si se considera un 
fluido incompresible (agua) dentro del dominio 
2 C Ru durante un intervalo de tiempo [0, 1), la 
dinámica del flujo en cada punto (x, y) en un 
instante específico t está determinado por las 
variables de estado, densidad de masa p(x, 1), 
el campo de velocidad u (x, 1) y su energía e (x, 
1). El planteamiento diferencial se deriva de la 
aplicación de los principios de conservación de 
masa, momento y energía sobre un volumen de 
control, características incluidas en las ecuacio- 
nes de Navier-Stokes(N-S), a partir del cual es 
deducida en su forma diferencial. Las ecuacio- 
nes de continuidad y cantidad de movimiento 
están expresadas por las ecuaciones (3) y (4), las 
cuales se generaliza para x, y, 2 (White, 1994; 
Cengel de Cimbala, 2012) 


Ecuación de continuidad: 


du, 00, a 


a ay 


6) 


Cantidad de movimiento en x: 
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Donde 2, 2, du 


ax ay 
métrica del volumen finito analizado (V = www) 
con respecto a los ejes x, y y 2; p es la densidad 


es la deformación volu- 


(kg m*); ves la viscosidad cinética (> Jevar- 


P 
do yes la viscosidad dinámica; t es el tiempo (5), 


y ges la aceleración de la gravedad. 

El objetivo del presente trabajo consiste en 
aportar información sobre el funcionamiento de 
un dispositivo Venturi bajo criterios de ecuacio- 
nes deducidas en forma diferencial e integral, 


partiendo de un dispositivo con caracter 
específicas y datos experimentales de laborato- 
rio reportados por Sotelo-Ávila (1979), así como 
comparar los resultados experimentales con los 
obtenidos por simulación usando DFC para de- 
terminación de parámetros de funcionamiento 
e instalación en sistemas de riego. 


¡cas 


Materiales y métodos 
Venturi en estudio 


El dispositivo descrito por Sotelo-Ávila (1979) 
corresponde al estudio reportado por Smetana 
(1957) y esquematizado en la figura 1. Tiene un 
diámetro en la zona inicial de D,= 0-15 m y un 
diámetro en la zona de estrangulamiento de 
D,=0.075 m. La transición entre estos diámetros 
en la zona de contracción es de forma cónica. 
La zona de expansión consta de una ampliación 
que va desde el final de la zona de contracción 
hasta alcanzar de nuevo el valor del diámetro 
D,. Según Sotelo-Ávila (1979), para su análisis 
se propuso un ángulo (0) de 25 a 7” como de- 
terminante de la longitud necesaria de la am- 
pliación de salida y con ello lograr normalizar 
el flujo. El dato en la deflexión del manómetro 
diferencial entre los diámetros D, y D, es de 
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A Distribución de presion 


PRES 


Ds E Arg 


Entrada 


Contracción: 


Expansión. 


Figura 1. Venturímetro en una tubería con Aujo a presión. 


Al: = 0.15 m de Hg, tal como se muestra en la 
figura 1; siendo el peso específico del mercurio 
ts, = 13 595 kg mv y la viscosidad cinemática del 
agua a 10C, y=0.013 cm's”. 

Las principales variables que se estudi 
este trabajo, por su influencia en el funciona- 


nen 


miento del aforador Venturi, son las corrientes 
cinéticas. En efecto, un diseño impropio de las 
zonas de entrada o salida se vería reflejado por 
una velocidad de flujo excesiva, la generación 
de zonas muertas o bien turbulencias. Según 
Webber (1971), el diámetro de la contracción 
D, debe ser lo suficientemente grande como 


para que la presión no caiga por debajo de dos 
metros de columna de agua con respecto a la 
presión absoluta, pues para esta presión hay 
tendencia a la formación de burbujas de aire, 
existiendo el peligro de cavitación en las pare- 
des del dispositivo. 

El Venturi mostrado en la figura 1 cuenta 
con un manómetro diferencial de mercurio, 


La observación experimental en laboratorio 
proporciona el valor de At y con empleo de 
la ecuación (5) (Sotelo-Ávila, 1979) es posible 
determinar el gasto (Q): 


65) 


Donde se denota por A, al área de la sección 
contraída (m?); g, la aceleración de la gravedad 
(m s?) y,,, como peso específico del mercurio 
(kg m”); y, como peso específico del agua (kg 
m); Ah, como la deflexión del manómetro de 
mercurio (m), y Cd como el coeficiente de des- 
carga (adimensional). 

En el cuadro 1 se indican las características 
geométricas del aforador reportado por Sotelo- 
Ávila (1979) y de acuerdo con Smetana (1957), 
las longitudes de entrada y salida son mayores 
que cuatro veces el valor de D,. 


Cuaino 1. Características geométricas del Venturi en estudio (Sotelo-Ávila, 1979). 


Dime) A ll E 
Radio | Ángulo 
2 4D, . 18 ly | grado | PP 
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Luego entonces se busca determinar el gasto 
a partir de los resultados experimentales para el 
aforador Venturi con las caracterí 


ticas ya seña- 
adas. Para ello, con la información reportada 
por Sotelo-Ávila (1979), se consideró como valor 
de velocidad inicial V, =1.533 m s'; como límite 
de estrangulamiento un valor de m1 = 0.25, con 
los cuales se mantenían valores del número de 
Reynolds mayores que 1 x 10”, Cabe resaltar que 
para la instalación de los aforadores Venturi en 
sistemas de riego, los parámetros de operación 
importantes que deben conocerse son: la pérdi- 
da de energía total generada en el dispositivo y 
la presión absoluta en la contracción. De acuer- 
do con Webber (1971), se recomienda realizar 
la medición de la línea de energía desde una 
distancia de 4D, antes de la contracción hasta 
4D, después de la ampliación, justo cuando el 
flujo se comporte como un fluido totalmente 
desarrollado. 


Dinámica de fluidos computacional (DFC) 


Tal y como fue mencionado, herramientas nu- 
méricas como la DEC han resultado eficaces en 
el diseño y simulación de flujos con base en una 
modelación a través de métodos aproximados 
de solución de las ecuaciones gobernantes de 
la mecánica de fluidos (ecuaciones (1), (3) y (4)) 
y esquemas numéricos de solución (volumen 
finito), con sus correspondientes algoritmos de 
procesamiento. Así, la DEC sigue una serie de 
pasos que pueden englobarse en tres grandes 
etapas: 


a) Preproceso. Consiste en la generación de la 
representación digital en tres dimensiones 


de la geometría del modelo físico del dis- 
positivo en estudio (figura 2a). Este paso 
es indispensable para definir el mallado de 
análisis (figura 2b), en el cual el dominio 
computacional es subdividido en un núme- 
ro finito de elementos prismáticos, donde 
las ecuaciones de conservación de energía, 
masa y cantidad de movimiento lineal serán 
aplicadas. 

b) Proceso. Una vez generada la geometría, 
se importa al procesador, donde se realiza 
la adjudicación de valores a las caracte- 
rísticas de frontera (figura 2c), entrada 
lado izquierdo y salida lado derecho; de 
igual manera, se establecen las hipótesis 
correspondientes respecto al Solver a uti- 
lizar, para conseguir la solución, basado 
en las ecuaciones de Reynolds Averaged 
Navier-Stokes Simulation y Presure-Based, 
indicando que se resuelve por separado 
la ecuación de Momentum para U, 0, 10, 
seguida de la ecuación de la energía; des- 
pués se solventa el modelo de turbulencia 
K-e, modelo más ampliamente utilizado e 
incorporado al software DFC los atributos 
del procedimiento, se indica en las cons 
deraciones (cuadro 2). 

e) Postproceso. Los resultados se muestran 
de manera gráfica y escalar (figura 2d). 
Cabe señalar que en el presente trabajo, 
la herramienta de DFC utilizada ha sido 
el programa comercial ANSYS Workbench 
v. 14.5 (Fluent, 1998), y cuenta con las in- 
terfaces necesarias para la generación de 
geometría del dispositivo (Design model), el 
mallado (Meshing) y el simulador numérico 
correspondiente (Fluent). 


Cuadro 2. Condiciones de inicio e hipótesis de cálculo. 


Hipótesis de simulación 


Soler Presure-Base 
Formulación Implícita 

Condición de tiempo Estacionario 
Modelo de turbulencia Ke de dos ecuaciones 


| Segregado 


Escenarios de simulación 


Velocidad del fluido 
Escenarios de simulación 


Constante (velocidad) 
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2b) 
lo ] 


Figura 2 Construcción del modelo computacional por etapas. 


Resultados y discusión 


Método de revisión para el Venturi probado en 
laboratorio 


De acuerdo con la metodología propuesta por 
Sotelo-Ávila (1979) y la información de labora- 
torio sobre la velocidad de inicio 1, = 1.533 
m s”, se determinó el coeficiente de descarga 
C, (figura 3). Este valor se obtuvo con base 
en la relación del estrangulamiento mi = 0.25 
(dato, cuadro 1) y del coeficiente de descarga 
C,= 1.009. La observación en laboratorio sobre 


la altura de columna de mercurio en la deflexión 
del manómetro diferencial reportado por Sotelo- 
Ávila (1979) es igual a Al: = 0.15 m, y con la 
ecuación (5) se determinó el gasto Q = 0.02698 
m' 5. Así, con la ecuación de conservación de 
masa, se determinó la velocidad en la sección de 
contracción V, = 


13 ms, Por último, con la 
D, 


ecuación del número de Reynolds R 
obtuvo un valor de Re = 3.54 x 10%, con lo cual 
se cumplió con la recomendación de ser mayor 
que un valor de Re=1 x 10%. 


Relación grado de estangulamiento-coeficiente 


PANA 


Cofciete de descarga 


000 olo odo 0% 0% 0 
¡Grado de estrangulamiento 
cocoa —= ido poto 


odo od 


Figura 3. Coeficiente C,con relación al grado de 
estrangulamiento 


La figura 3 sustituye a las reportadas por 
Smetana (1957). Esta gráfica tiene por objeto 
mostrar el modelo matemático que representa 
la relación existente entre el grado de estrangu- 
lamiento (m) y el coeficiente de descarga (C), 
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destacando que el modelo encontrado cuenta 
con un coeficiente de determinación R*=0.9972. 

Con base en el procedimiento experimental 
convencional descrito, se observa que cada vez 
que se efectúe algún cambio en la velocidad 
inicial en el punto de entrada al Venturi, es ne- 
cesario medir el cambio correspondiente en la 
altura de la columna del manómetro diferencial 
de mercurio y con ello poder calcular el gasto en 
las nuevas condiciones del flujo. Es decir, el gas- 
to está relacionado con el valor Alt de la columna 
del manómetro de mercurio, lo que conduce a 
una sola ecuación, pero con dos incógnitas. Por 
lo anterior, con este procedimiento experimental 
es necesario conocer la velocidad inicial (en 
la entrada al dispositivo), el valor de Al: del 
'manómetro diferencial y el valor de la relación 
de estrangulamiento mm, para luego obtener 
de forma gráfica el valor correspondiente del 
coeficiente de descarga C, y por último calcular 
el valor del gasto. El principal inconveniente del 
procedimiento experimental descrito consiste 
en que cualquier cambio en las condiciones del 
flujo requiere del montaje experimental y de los 
datos observados. 


Método de revisión para el Venturi con DFC 


La figura 4 muestra los resultados obtenidos por 
simulación con DFC para la magnitud de los 
vectores de velocidad del flujo para cada punto 
(identificados por círculos) en las secciones 
de las zonas de entrada al Venturi en el punto 
D,, de contracción, de ampliación y de salida 
del dispositivo. Las magnitudes de velocidad 
están representadas por una escala gráfica de 
colores, En la figura 5 se muestran los valores 
de velocidad y presión puntuales del flujo y se 
relacionan con su ubicación en la sección del 
tubo en la zona de entrada; es decir, un valor 
nulo de velocidad en las paredes del dispositivo 
hasta un valor máximo de V =1.8 m s' en el 
centroide de la sección transversal del tubo y 
valores de presión de P =78.7 cm de Hg en la 
pared del tubo y de P =78.4 cm de Hg en el 
centroide de la sección transversal del tubo. Con 
esta información se determinó la velocidad y 
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no a y vay a a 


Figura 4. Distribución de velocidades. 


presión medias del flujo en el punto de entrada 
con ayuda de las ecuaciones (6) y (7): 


y VA VA +. +VA, 6 


ñ A 


Donde V,, = velocidad media del flujo en 
la sección; V, = velocidad puntual en el círcu- 
lo 1 hasta 11, A, = área en el círculo 1 hasta 11; 
A, = área total de la sección transversal del 


tubo: 


(7) 


DondeP,, 'n media del flujo en la sec- 
ción; P, = presión puntual en el círculo 1 hasta 
1; A, = área en el círculo desde 1 hasta 1 A, 
= área total de la 

Este procedimiento se repite para las seccio- 
nes correspondientes a las zonas de contracción 
y salida del dispositivo. Los valores de veloci- 
dad promedio de flujo presentes en las zonas 
del dispositivo estudiado se presentan en el 
cuadro 3. 

Al comparar los valores reportados experi- 
mentalmente por Sotelo-Ávila (1979) y los obte- 
nidos mediante DFC (cuadro 3), con base en el 
cálculo del error relativo absoluto |e, | (ecuación 
(8)), se observa que son prácticamente iguales. 
Es decir, una velocidad media en la zona de 


ción transversal del tubo. 
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Relación diámetro-velocidad-presión 
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020 se 78350 E 
0000 1 78300 
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Piámatro del tubo 
—otocidnd === Presión 
Figura 5. Peres de velocidades y presión (zona de entrada) 
Cuadro. Resultados e a simulación con DFC del Venturi experimental 
ENS PETT Zona de salda 
presión media | vetociónd Presión vetocitad | Femina — [ietociad 
E = media E media 
O O O O IA 
1006] 7 [ar [as [ww] os | ss [ws] 7 [ros [13 


entrada (153 ms .013), un gasto (0.027 
ms”, por conservación de masa); | e, | =0.013) y 
una presión en el manómetro diferencial, 15 cm. 
de Hg (diferencia de presión en cm de Hg entre 
la zona de entrada hasta la zona de contracción; 


alor observado - valor simulado 


(8) 


valor observado 


Una vez validado el modelo construido por 
DEC, se simularon ocho escenarios adicionales, 
con velocidades medias desde 0.028 hasta 2. 264 
ms, correspondientes a una variación de gasto 
desde 0.005 hasta 0.040 m' 5”. En el cuadro 4 se 
reportan los resultados de los ocho escenarios 
simulados (E2 a E9), donde la presión se expresa 


en pascales, cm de Hg y en kg cm; así como 
la velocidad media en las zonas de entrada, 
contracción y salida, en m s”. 

Con base en el análisis de los resultados 
mulaciones para el rango de 
gasto de operación del Venturi (Q = [0.005 m?s; 
Q=0.040 mv 5"), y reportados en el cuadro 4, 
se ha construido un modelo matemático de tipo 
regresivo entre las variables gasto del Venturi 
(Q en m' s*) y pérdida de carga en la zona de 
contracción del dispositivo (Alt en cm de Hg). 
Este modelo está representado por una ecuación 
de segundo grado y cuenta con un coeficiente 
de determinación R?= 1 (figura 6). Lo anterior 
es totalmente congruente con la génesis del 
fenómeno hidráulico local en la sección de con- 
tracción, pues la caída de presión se debe a un 
cambio drástico de velocidad del flujo y donde 


obtenidos en las 


nom. 
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Cuadro 4. Resultados de las simulaciones con DEC para los acho escenarios propuestos, 


ETA AA E 
Frasón —— | Velocidad Psión Velocidad Prosión Velecidad 
ión | rasctos| [ago me | pjs [dm ajo | gone] ed 
E [orar | 761 [100 | 0285 | moss [756 | 200 | 100 [100005 | 760 [100 
E [ro | 76 [100 0577 | 0930 05 [300 | 2062 [0030 | 760 108 
5s [rozas | 768 100 0506 | o6m1 [720 [099 | 22 [00322] 760 108087 
5 | ma [15 | 130 | 0927 [05 [095 | 00 [0037] 760 100 im 
ss [rosa | 782 [106] 1020 | 00008 [667 [091 | 5160 [0132] 760 [101 
E [53m | mo 108 | 175 | somo [01 000 | ca [00m | 760 [10010 
Es [106700] 401 109 | 1093 | r030 [573 07 sora | 760 100208 
5 Ja] m2 [10] 220 | 0005 [515 [om] ss [00309] 760 [100] 22 
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Pérdida de carga en em de mercurio 


o pl cin — 


Figura 6. Determinación de las pérdidas en la sección 


contracción a partir del valor del gasto. 


las pérdidas locales por fricción también están 
en función de la velocidad del flujo. Con este 


modelo se reduce la incertidumbre, al eliminar 


el requerimiento de selección del coeficiente de 
descarga C, en el cálculo del gasto ante cual- 


q 
modelado. Es decir, con ello se ha simplificado 


» otra velocidad inicial dentro del rango 


el problema a una ecuación con una sola incóg- 
nita. Por ejemplo, usando la gráfica de la figura 
6 y considerando el diámetro del tubo del caso 
de estudio (D,=0.15 m), si se tiene una lectura 
en el manómetro diferencial de presión de 
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Al: =0.15 m de Hg, se obtiene un gasto Q =0.027 
* s!, que corresponde a una velocidad de 1. 


Por otro lado, con la información proporcio- 
nada por las simulaciones con DFC, es posible 
construir un modelo matemático adicional que 
facilite la estimación de la pérdida total de carga 
ocurrida en el aforador Venturi para el rango de 
gastos estudiados. En efecto, en forma equiva- 
lente al caso anterior, se considera el valor de 
caída de presión entre la sección inicial y final 
del dispositivo, relacionando este diferencial 
de presión con el gasto de entrada al Venturi. 
El modelo encontrado presenta las pérdidas de 
presión (Ah) en cm de columna de H,O versus 
el gasto (Q en m* 5") mostrado por la figura 7, 
y cuenta con un coeficiente de determinación 
(R) superior a 0.99. Cabe resaltar que este tipo 
de modelo no es posible construirlo cuando 
sólo se cuenta con estudios convencionales de 
experimentación en laboratorio. 

Así, por ejemplo, utilizando el modelo de 
la figura 7, para un gasto máximo de 0.040 
m' s? se obtiene una pérdida de carga total en 
la longitud de Venturi de 0.075 kg cm”. Esta 
pérdida total de carga está considerada entre 
la distancia formada por + 4D,, a partir de la 
sección de contracción y donde D, es la sección 
de entrada al dispositivo. 


Iniguez-Covarrubias, Comparación de resultados experimentales de un Ventura con simulación de dinámica de fluidos computacional 


Relación gasto-pérdidas totales 
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Figura 7. Determinación de las pérdidas totales ocurridas 
enel dispositivo Venturi entre las secciones inicia y final, 


Finalmente se resaltan otros tres resultados 
importantes que sólo pueden derivarse con base 
en la simulación efectuada con DFC: a) la carga 
hidráulica de funcionamiento en la contracción 
del Venturi deberá respetar una presión mínima 
disponible de 0.4 kg em” para el gasto máximo 
de 0.040 m' s* y con ello satisfacer la condición 
sugerida por Webber (1971) para evitar la gene- 
ración del fenómeno transitorio de cavitación; 
b) por la condición anterior, el gasto límite 
superior para las dimensiones y condiciones 
de trabajo es precisamente el gasto de 0.040 nv 
1; y e) la pérdida de carga total para la longitud 
de Venturi es de 0.075 kg cm para Qs, =0.40 


Conclusiones 


En este trabajo se modeló, con ayuda de la 
dinámica de fluidos computacional, el funcio- 
namiento de un dispositivo de medición de 
Aujo tipo Venturi. La información experimental 
obtenida en laboratorio y reportada por Sotelo- 
Ávila (1979) facilitó la construcción del modelo 
matemático y su validación, obteniéndose 
una muy buena congruencia entre los valores 
experimentales y los simulados de velocidad 
del flujo, gasto y caída de presión en el manó- 
metro diferencial. La modelación con DFC fue 
realizada con el empleo del programa Fluent, 
el cual utiliza el método del volumen finito y 


se encuentra inserto en el software comercial 
ANSYS Workbench u. 14,5. 

Con base en lo aquí analizado es posible 
enumerar las conclusiones siguientes: 

1. Existe una congruencia altamente satis- 
factoria entre los resultados medidos en forma 
experimental en laboratorio y los encontrados 
por simulación numérica con DEC. 

2. Los resultados obtenidos de los esce- 
narios hidráulicos modelados con DEC han 
permitido comprender con mayor detalle el 
funcionamiento del Venturi y la determinación 
de modelos matemáticos de tipo regresivo sobre 
variabilidad de sus características de operación. 
Los modelos encontrados a partir de los resul- 
tados con DFC permiten, con un alto nivel de 
precisión, la determinación de condiciones de 
operación para situaciones no experimentadas 
ni simuladas. 

3. En particular, con el procedimiento aquí 
propuesto, ha sido posible modelar, a través 
de una ecuación de segundo grado, la relación 
gasto versus pérdida de carga en la sección de 
contracción del Venturi y con ello eliminar la 
incertidumbre en la elección del coeficiente 
de descarga C,. Igualmente fue modelada la 
relación gasto versus pérdida de carga total en 
el dispositivo Venturi, característica difícil de 
modelar bajo un contexto tradicional de expe- 
rimentación en laboratorio. 

4. La modelación con DFC puede ser usada 
para el diseño y la definición de dimensiones 
óptimas de un Venturi, considerando diferentes 
tamaños, formas, condiciones de funcionamien- 
to y materiales, además de poder determinar las 
condiciones límite que ocasionarían la presencia 
del fenómeno transitorio de cavitación en las 
paredes del Venturi 

Por último, es importante resaltar que la 
modelación con DFC y la experimentación 
en laboratorio son métodos complementarios 
de análisis de fenómenos hidráulicos, y su 
combinación es la mejor estrategia para el 
encuentro de una solución técnica, práctica, 
económica y confiable, en particular en el caso 
de flujómetros tipo Venturi y su instalación en 
sistemas de riego. 


nom. 


isSNoLgr306 


Ihiguez-Covarrubias, Comparación de resultados experimentales de un Venturi con simulación de dinámica de fuidos computacional 


Referencias 


Anderson. D. Jr (1995), Computational Fluid Dynamics. The 
Basics with Applications (328 pp. New York: MeGraw- 
E 

ASME (1983), Fluid Metrs, this Theory and Application (Gh 
edition). In H.S. Bean (Ed.) New York: American Socioty 
of Mechanical Engineers, Research Commitoe on Fluid 
Meters 

Gengel, Y. A, Se Cimbala,J. M. (2012), Mena de fividos- 
fundamentos y aplicaciones (978 pp.) (2a edición). México, 
DF: McGraw Hill 

Chow, Y. T. (1959). Open Channel Hydraulics. New York 
MeGraw-Hil 

Davies, J. A, £ Stewart, M. (2002) Prctng Globe Control 
Valoe Performance. Part1. CFD, Modeling. ASME Journal of 
Ftuids Engineering, 12465), 772-777. 

Flores-Velazquez, J, Mejia, E, Rojano, A, €: Montero, 
J. (011). Análisis del clima en un invemadero con 
ventilación mecánica, Agrociencia,45(5), 545-560. 

Fluent (1998). Fluent, 5. Sheffcld, UK. Fluent, Europe Lt 

Gharahjeh,S., Aydin, L, de Altan-Sakarya, A. B. (2015), Weir 
Velocity Formulation for Sharp-Crested Rectangular 
Wir, Fox Measurement and Tostrumentation, 41, 50-56. 

1OS (1991). Mensurement of Flow by Means of Pressure 
Diferential Devices. Geneva: International Organization 
forStandardization ISO 51674. 

Levy, E (1957). Mexúnica delos fluidos (266 pp). México, DF: 
Instituto de ingeniería, UNAM. 

Manzano, ], de Palau, G. (2005). Hudraulic Modeling o Venturi 
únjector y MemsofCFD. Tampa, USA: ASAE International 
Meeting. 

Naghavi, B, Esmaili, K, Yazdi J, 8e Vahid, E K. (2011). 
An Experimental and Numerical Study on Hydraulic 
Characteristics and Theoretcal Equations of Circular 
Wir, Canadian Journal of Cil Engineering, 38(12), 1327. 
1334 

Palau, G, Arviza,], de Frankel,S, (2004). Tre-Dimensional 
Control Valve with Complex Geometry: CFD Modeling and 
Experimental Validation. 34TH AAA. Fluid Dynamics 
Conference and Exhibit. IN AIAAA-2001- Portland, 
USA. 

Ranald, V. G. (1978), Mecónico de fuidos e Hidráulica (2a 
edición), México, DF: MeGraw-Hil 

Smetana J. (1957). Hidrdulka. Praga: Ceskoslovenka 
Akademie, VED. Referencia de Sotelo-Ávila, G. (1995) 
Hidráulica general (277 pp.) México, DF Limusa. 

Sotelo-Ávila, G. (1979). Hidnúulica general (277 pp.). México, 
DF: Limusa. 


ISSN DIS73%6 


SARH (1973). Proyecto de Zonas de Riego (567 pp). México, 
DE Dirección de Proyectos de Irigación, Departamento 
de Canales. 

Temecpattanapongsa, S, Merkdey, G. P, Barfuss, S.L, e 
Smith, B. (2014). Generic Unified Rating for Cutroat 
Fltmes,Ivigaton Science, 32) 29-40. 

USDI (1979). Water Measurement Manual (3rd. edition) 
Denver: US Department of Interior Burcaw of 
Reclamation, US Govern Printing Office. 

Webber, NB, (1971). Fluid Mechanie for Civil Enginsers 
pp). London: Chapman de Hall 

White E. (1994). Fluid Methamios (736 pp.) 3rd edition). New 
Yorke McGraw Hill Inc 

Wu, W, Rodi, W, £e Wenka, T. (2000). 3D Numerical 
Modeling af Flow and Sediment Transport in Open 
Channels. ASCE Journal of Hydrauic Engineering, 12600) 
45 

Zheng G. (2000). Applications of CED Tools to Design and 
Development af Pumps. Turbomchiery, 2811) 649-656. 


30 


Dirección institucional de los autores 


Dr. Mauro Iñiguez Coarrubias 
Dr. Jorge Flores Velazquez 

Dr. Waldo Ojeda Bustamante 
Dr. Roberto Mercado Escalante 


Instituto Mexicano de Tecnología del Agua: 
Paseo Cuauhnáhuac 8532, Col. Progreso. 
162550 Jiutepec Morelos, México 

Teléfono: +52 (777) 329 3600 

Fax: +52(777) 319 4220 
micitlalocimta.mx 

jorge MoresGtlalocimtamx 
wojedaWHalocimta.mx 
rmercadofHaloc.imta.mx 


Dr. Carlos Díaz-Delgado 


Universidad Autónoma del Estado de México 
Centro Interamericano de Recursos del Agua 
Facultad de Ingeniería 

Carretera Toluca Atlacomulco kim 145, Unidad San 
Cayetano 

50200 Toluca, Estado de México, México. 

Teléfono: +52 (722) 2965 550 

ediazdiluaemexemx 


